
X=-1+2.*rand(1,10);
Y=-1+2.*rand(1,10);

szX=size(X,2);
testpos=[];

for i=1:szX
    for j=1:szX
        if (X(i)*X(i)+Y(j)*Y(j))<1
            testpos=[testpos;[X(i),Y(j)]];
        end
    end
end
Hm=[];
for i=1:size(testpos,1)
    Hx=0.0;
    for j=1:T
        if(features(j,1).*testpos(i,1) +features(j,2)>testpos(i,2))
            Hx=Hx+features(j,5);
        else
            Hx=Hx-features(j,5);
        end
    end
    Hm=[Hm,Hx];
end
Hm